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Abstract. First we prove a general spectral theorem for the linear Navier- 
Stokes (NS) operator in both 2D and 3D. The spectral theorem says that the 
spectrum consists of only eigenvalues which lie in a parabolic region, and the 
cigcnfunctions (and higher order eigenfunctions) form a complete basis in H e 
{I = 0, 1, 2, • • • ). Then we prove the existence of invariant manifolds. We are 
also interested in a more challenging problem, i.e. studying the zero-viscosity 
limits [y — > + ) of the invariant manifolds. Under an assumption, we can 
show that the sizes of the unstable manifold and the center-stable manifold 
of a steady state are 0(\/v), while the sizes of the stable manifold, the center 
manifold, and the center-unstable manifold are 0{v), as v — > + . Finally, we 
study three examples. The first example is defined on a rectangular periodic 
domain, and has only one unstable eigenvalue which is real. A complete esti- 
mate on this eigenvalue is obtained. Existence of an ID unstable manifold and 
a codim 1 stable manifold is proved without any assumption. For the other 
two examples, partial estimates on the eigenvalues are obtained. 
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1. Introduction 

Navier-Stokes equations define an infinite dimensional dynamical system. Tur- 
bulence can be represented through a chaotic solution in the infinite dimensional 
phase spaces. For example, one can take the Sobolev spaces H e {I — 0, 1, 2, • • • ) to 
be the phase spaces. A basic dynamical system question is: Does the phase space 
foliate under the NS flow ? That is, are there invariant manifolds under the NS flow 
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? The answer is yes. It seems that under the Euler flow, which can be obtained 
from the NS flow by setting viscosity to zero, the phase space does not foliate. To 
begin a dynamical system study, one chooses a common steady state of the NS and 
its corresponding Euler equation, defined on a periodic spatial domain. At this 
steady state, the spectra of the linear NS and the linear Euler are dramatically 
different. The linear NS has the form: i/A+ relatively compact terms, therefore, 
it has only a point spectrum via Weyl's essential spectrum theorem. On the other 
hand, the linear Euler has a nontrivial essential spectrum. For 2D linear Euler, if 
the steady state as a vector field has a saddle, then the essential spectrum in the 
Sobolev space H l (£ = 0, 1, 2, • • • ) is a vertical band symmetric with respect to the 
imaginary axis |25J. The width of the band (= £A for some A > 0) increases with 
£, representing the fact of cascade and inverse cascade. A typical example of such 
steady states is the cat's eye. If the steady state has no saddle, but has at least two 
centers, then the essential spectrum in the Sobolev space H is the imaginary axis 
|25j . A typical example of such steady states is the shear |16| |9]. For the shear, 
a sharp upper bound on the number of eigenvalues was also obtained [H]. For 3D 
linear Euler, the spectrum is an open problem. A nice example of the steady states 
is the ABC flow vector field. 

An interesting open problem is to understand the connection between the spec- 
tra of linear NS and linear Euler as v — > + . Besides the progresses mentioned above 
in understanding the spectra of linear Euler, progresses were also made in better 
understanding of the spectra of linear NS, although they were in a different setting. 
For both 2D and 3D NS posed in L 2 (in terms of weak solutions) on a bounded 
domain with non-slip boundary condition, Prodi |22| proved that the eigenvalues 
of the linear NS at a steady state lie inside a parabolic region X r < a — bXf where 
a > and b > are some constants, and A = A r + i\ is the spectral parameter. 
Moreover, David Sattinger 24 proved that the linear NS has infinitely many eigen- 
values lying in the parabolic region, each of finite multiplicity, which can cluster 
only at infinity, and the corresponding eigenfunctions (and higher order eigenfunc- 
tions) form a complete basis in L? . Similar results are proved in the current article 
for a periodic domain and mild solutions in H {I = 0, 1, 2, • • • ). The proof involves 
some deep result from Hilbert-Schmidt operators. Together with some deep result 
from sectorial operator, we prove the existence of invariant manifolds. 

In terms of well-posedness of Navier-Stokes equations (NS), zero- viscosity limit 
problem has been extensively studied |7] |26j . Here we are also interested in the 
zero-viscosity limit (y — > + ) problem for invariant manifolds. It is proved in this 
article that under an assumption, their sizes are of order 0{yfv) or 0(y) as v — > + . 
Existence of invariant manifolds for Euler equations is still open. I tend to believe 
that invariant manifolds for Euler equations do not exist, at least not in the usual 
sense. An intuitive argument is as follows: Kato's technique |2j for proving the 
existence of a mild solution to Euler equations amounts to writing the solution as 
u(t) = e~ Jo tl ( T )' v<ir u(0), where u is the velocity, and the Leray projection and 
the external force are not included, which do not affect the argument. Notice that 
the u(t) in the exponent depends upon u(0), therefore, one should only expect the 
evolution operator of Euler equations to be at best C° in u(0), just as in t. On the 
other hand, the existence of an invariant manifold is really an at least C 1 in u(0) 
phenomenon. One can obtain expressions of invariant manifolds around a shear for 
a Galerkin truncation of the Euler equation |10| . which exhibit a lip structure. 



INVARIANT MANIFOLDS FOR NAVIER-STOKES EQUATIONS 



3 



For concrete examples, e.g. single mode steady states, a continued fraction 
technique can be designed for calculating the eigenvalues |16| . Detailed calculations 
were conducted by Vincent Liu |13| |14| . Liu obtained rather detailed information 
on unstable eigenvalues. In some case, Liu showed that as v — > + , the unstable 
eigenvalue does converge to that of linear Eulcr. In this article, besides the example 
of Liu, we will study two more examples. For the first example, a complete eigen- 
value information is obtained, and the existence of unstable and stable manifolds is 
proved without any assumption. For the second example, only partial information 
on the eigenvalues is obtained, as in the example of Liu. 

The article is organized as follows: In section 2, we will present the formulation 
of the problems. Section 3 is on invariant manifolds. Examples are studied in 
Section 4. 

2. Formulation of the Problems 

We will study the following form of 2D Navier-Stokes equation (NS) 

(2.1) <9 t ft + {^,ft} =i/[Aft + /(a:)] , 

where ft is the vorticity which is a real scalar-valued function of three variables t 
and x = (x\, x%), the bracket { , } is defined as 

{f,9} = (d x J)(d X2 g)-(d x J)(d xl g) , 

where $ is the stream function given by, 

ui = -9 Xa * , u 2 = d Xl ^ , 

the relation between vorticity ft and stream function \& is, 

ft = d Xl u 2 - d X2 u x = A* , 

and v is the viscosity, A is the 2D Laplacian , and f(x) is the body force. We will 
pose periodic boundary condition 

ft(i, xi + 27r, x 2 ) — ft(t, xi,x 2 ) = ft(t, xi,x 2 + 2ir), 

i.e. the 2D NS is defined on the 2-torus T 2 . We require that if? and / have mean 
zero 

/ Vdx = I fdx = . 

Of course ft always has mean zero. In this case, $ = A _1 ft. 

Setting v = in the 2D NS (|2.1|l . one gets the corresponding 2D Euler equation. 
We will study the common steady states of 2D NS and 2D Euler given by 

(2.2) {*,ft} = 0, Aft + /(a;) = 0. 

Our three prime examples are: (1). The shear ft = cos 22 defined on the rectangular 
periodic domain [0, 2ir/a] x [0, 27r] where 1/2 < a < 1. (2). The shear ft = 
cos(a:i + x 2 ) where fix) = 2cos(iei + x 2 ). (3). The cat's eye ft = cosxi + acosa^ 
where a is a constant and f(x) = cosa;i + a cos x 2 . 

We will study the following form of 3D Navier-Stokes equation 

(2.3) <9 t ft + (u ■ V)ft - (ft • V)it = ^[Aft + f(x)] , 

where u — [u\, 112,113) is the velocity, ft = (fti,ft2,fta) is the vorticity, V = 
(d Xl , d X2 , d X3 ), ft = V x u, V ■ u = 0, v is the viscosity, A is the 3D Laplacian, 
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and f(x) = (/(xi), /(X2), /(X3)) is the body force. We also pose periodic bound- 
ary condition of period (2tt, 2tt, 2tt), i.e. the 3D NS is defined on the 3-torus T 3 . 
We require that u, f and f2 all have mean zero. In this case, u can be uniquely 
determined from fl by Fourier transform: 

Mk) = i\k\- 2 [k 2 n 3 (k) - k 3 n 2 (k)} , 
u 2 {k) = i\k\- 2 [k 5 h x {k)-k x h s {k)}, 

u 3 (k) = i^-^k^k) ~ k^k)] . 

Setting v = in the 3D NS (|2.3[1 . one gets the corresponding 3D Euler equation. 
We will study the common steady states of 3D NS and 3D Euler given by 

(2.4) (tt- V)0- (O- V)u = , AQ + f(x)=0. 

Our prime example is: The ABC flow 

ui = A sin £3 + Ccosx2 , u 2 = Bsmxi + Acosx 3 , u 3 = C smx 2 + B cosxi . 
Here SI = u = f. 

3. Invariant Manifolds 

Centering around a C°° steady state f2*, it*, 'J'* given by l|2.2|) or l|2.4|) . one 
can rewrite the NS as 

(3.1) d t fl = LVL + N(Q) , 

where for 2D 

LVt = vAfl ~ , 
N(Q) = -{^,0}, 

and for 3D 

LVL = uAQ - (u* ■ V)0 - (u ■ v)n* 
+ • V)it + (O • V)u* , 
iV(f2) = -(tt- V)Q+ (n- V)u . 

Theorem 3.1. T/ie linear Euler operator L is a closed operator. Its spectrum 
consists of infinitely many discrete eigenvalues, each of finite multiplicity, which 
can cluster only at infinity. All the eigenvalues lie inside the parabolic region 

X r < a — b\f 

where a > and b > are some constants, and A = A r + i\ is the spectral 
parameter. The corresponding eigenf unctions (and higher order eigenf unctions) 
form a complete basis in H (£ = 0, 1, 2, • • ■ ). 

Proof. The proof for L being a closed operator is trivial. We can write L as 

L = uA + A 

where A is relatively compact with respect to i/A. By Weyl's essential spectrum 
theorem |23j . L and i/A have the same essential spectrum. Since the essential 
spectrum of vA is empty, so is L. Therefore, L has only eigenvalues. It is easy to 
see that for large enough A > 0, (L — A) -1 is compact, thus each eigenvalue can 
be of only finite multiplicity and all the eigenvalues can cluster only at infinity [2] . 
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Next we show that large enough A > is not an eigenvalue. Assume the contrary, 
then there exists nonzero ft <G H t+2 such that 

\n - vAn = An . 

Applying (— A) £ / 2 , one gets 

X(-A) e / 2 n - vA(-A) e / 2 n = {-Af' 2 {AVt) . 
Multiplying by (— A) e / 2 Q, one obtains 



{-Af 2 n 



+ v 



(-A) £ / 2 d • \(-A)(-Ay/ 2 n 



(-A) £ / 2 d • \(-A) e / 2 (An) 



Thus 



(_ A )(*+1)/2 Q 



{-Af/ 2 n + 

= J [(-A) e / 2 n\ ■ [(-A) e / 2 (An) 
c J \(-A)^ 2 n\ \(-Af +1 ^ 2 n 

c[e 2 ||(-A)( £+1 )/ 2 0| 2 + 1 j\{-Af 2 ^ . 



c 



< 



Let e 2 = ^ and A = + 1 , we have 



J Y-Af' 2 ^ J \(-Af +1 ^ 2 n 



< 



which is a contradiction to f2 being an eigenfunction. Thus A = ^ — h 1 is in the 
resolvent set. Next let A be a complex number in the resolvent set. Naturally, 
we need to consider as a complex- valued function. For any h € H e , there is a 
Q e H t+2 such that 

-vAVt + \fl-Afl = h. 

Applying (— A) £ / 2 , one gets 

-isA(-A) e / 2 n + \(-Ay/ 2 fl - {-Af' 2 {Aa) = {-A) l ' 2 h . 
Multiplying by (-A) £ / 2 fj, one obtains 



(-A)(-A// 2 • (-A) £ / 2 



+ A 



{-A) l ' 2 n 



'-AY' 2 n 



{-Af 2 {An)] ■ \{-Af' 2 n\ + \{-Af 2 h] ■ \{-Af' 2 n 



Thus 



(3.2) 



( _ A) (^+D/2 



+ A 



{-A) l ' 2 n 



-I 



(-A)^ 2 (An) ■ {-AY' 2 n 



+ 



I 



{-A) l ' 2 h ■ (-A// 2 fi 
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Taking the real part, one gets 



v\\nWi +1 + \ r \\n\\ 2 e < c(e 2 \\n\\ 2 e+1 + ±\\n\\ 2 e 



+ ^(IWIl + IMU 2 ) 



Let e — 5jy, one has 



(3.3) oHI^+i + VllOIII^-H/illl + dllOII 



Taking the imaginary part of (|3.2I) , one gets 

A*||fi||!<C 2 ||fi|U +1 ||0|U + \\h\\ t \\Sl 



Xi\\n\U< c 2 \\n\\ e+1 + \\h\\ 



xm\ 2 i<c 3 Ml + i + c4h\ 12 



i > 



that is 

Then we have 
that is 

(3-4) ^llfill^^llMI^ + gll^ 2 

Adding (|3.3|l and (|3.4jl . we have 

V \2 \ IIOII2 - 1 A , ^4 



Ar + ^-A^ - dj ||n||i < - I 1+ ) ||fc|| 



Let a = C\ , 6 = and C 5 = i ( 1 



_ _ I/C4 

2C 3 ' ^ 2 ^ C 3 

(A r -a + 6A?)||0||| < C B ||/i||? 



Thus if A r - a + bXf > 0, 



(3.5) A) _1 || < 



a 



1/2 
5 



(A r -a + 6A2)V2 



(L — A) -1 is holomorphic in A and can be continued analytically as long as \\(L — 
A) _1 | is bounded. As shown above, (L — A) -1 exists at large A > 0. Then from 
(|3.5I) . (L — A) -1 exists everywhere outside the parabolic region A r < a — bX 2 . Thus 
the eigenvalues of L lie inside this region. To prove the rest of the theorem, we 
need a known lemma [Corollary 31, pp.1042, [2]]. 

Lemma 3.2. Let T be a densely defined unbounded linear operator in a Hilbert 
space H, with the property that for some Ao in the resolvent set of T , (T — Ao) _1 
is a Hilbert- Schmidt operator. Let 71, • • • ,75 be non-overlapping differentiate arcs 
having limiting directions at infinity, and such that no adjacent pair of arcs form 
an angle as great as ir/2 at infinity. Suppose that the resolvent (T — A) -1 satisfies 
an inequality \\(T — A) _1 || = 0(|A| _Ar ) as A — > 00 along each arc jj for some 
natural number N . Then the eigenf unctions (and higher order eigenfunctions) form 
a complete basis in H. 
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From lj3.5|l . we see that along any ray A r = c\i (c finite). 

\\{L-\)- 1 \\=0{\\\- 1 ) as A^oo. 

Next we show that (L — A) -1 is a Hilbert-Schmidt operator for large enough A > 0. 
Notice that 

(L - A)" 1 = [I + {vA - A) _1 A] _1 (i/A - A)- 1 . 



Since 



2( 



-y 



1 Ifcl 2 

\k\ 2 + \/v\k\ 2 + \/v 

<^EW 2( '~W, 



lAI^-^in* 



2 



we have 

It is easy to see that 
Thus 

Choose Ao = then 



II^A-A)- 1 ^!, < -±=||n||<_i 
V i/A 



\A£l\t-\ < c\\n\\ 



\\(pA- X)- 1 AQ\\ e < -^=||fi|| 
vVA 



I^A-Ao)- 1 ^!! < ^ 



We have [7 + (z/A — Ao) 1 A\ 1 being a bounded operator. Recall that a bounded 
linear operator T is a Hilbert-Schmidt operator if 

is finite where {e^} is a complete orthonormal basis [3]. By using the Fourier 
basis for H e , one sees that (i/A — Ao) -1 is a Hilbert-Schmidt operator. A product 
of a bounded linear operator with a Hilbert-Schmidt operator is another Hilbert- 
Schmidt operator. Thus (L — Ao) -1 is a Hilbert-Schmidt operator. By Lemma 
13.21 the corresponding cigcnfunctions (and higher order eigenfunctions) of L form 
a complete basis in H l (I = 0, 1, 2, • • • ). □ 

We will study the case as described by the following set-up. 

• Set-Up. Let there be m unstable eigenvalues A" ■ • • A^, n center eigen- 
values Af • • • A° , and of course the rest infinitely many stable eigenvalues 
Af A| • • • . Let the real parts of A" and Af have the smallest absolute values 
among the unstable and stable eigenvalues, respectively. 

Remark 3.3. For specific examples, often there is no center eigenvalue, here 
we put them in for generality. We need 77 (£ = 0, 1, 2, • • • ) to be a Banach algebra 
PQ. Thus for 2D and 3D, £ > 2. Our argument requires that we work with H e 
where £ > 3 for both 2D and 3D. 
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We have the splitting 

where iff (z — u, c, s) are spanned by the corresponding cigenfunctions. Projections 
onto Hi are denoted by P z (z = u, c, s). Denote by 

L Z =P Z L, N Z =P Z N, Sl z = P z n, (z = u,c,s). 

Then (|3.1[1 can be rewritten as 

(3.6) d t n z = L z n z +N z (n) , (z = u,c,s). 

Since in this case the spectral mapping theorem is trivially true, one has the tri- 
chotomy: 

(3.7) ||e tL "|| < Ce at , t < ; 

(3.8) ||e tLC || < C\t\ ni , teR; 

(3.9) ||e tLS || < Ce~ pt , t > ; 

where a = Re{A"} - e, f3 = -Re{Af} — e, m < n, and < e < -Re{Af}. Our 
argument depends upon the following two crucial facts: 

(1) One has a stronger inequality than l|3.9[l [[5], Theorem 1.5.4] due to the 
fact that L s is sectorial, 

(3.10) ||e tL "r! s ||, + i < Cm-^e-^We , t > . 

(2) Due to finite dimensionality, one has 

(3.11) m z h+i <C\\n z \\ t , (z = u,c). 

Remark 3.4. For (|3.11|) . take z — u for example, since fi" • • ■ belong to H l 
for any t, we can take 



(3.12) C= max \ / \\^\\ t 

l<j<m I 

Theorem 3.5. In a neighborhood of the fixed point fi» in the Sobolev space 
H e (T d ) (£ > 3, d = 2, 3), there exist a m- dimensional C°° unstable manifold, a 
n- dimensional C°° center manifold, a (m + n)-codimensional C°° stable manifold, 
a (m + n)- dimensional C°° center-unstable manifold, and a m-codimensional C°° 
center-stable manifold. 

Proof. We will take the center-stable manifold as the example. The argu- 
ments for the others are similar. Apply the method of variation of parameters to 
(|3.6(l . one gets the integral equations 

(3.13) n z (t) = e^-^Q'ito) + [ e {t - T ^N z (n)dr , (z = u, c, a) . 

J t 

To prove the existence of a center-stable manifold, we start with the Banach space 
B^ = \tt(t) ||0|| T ^ = supe-T* V \\n z (t)\\ e < +oo , (]a < 7 < a) \ . 

{ *>° z^c,s 4 J 

For such f2(i) in B 1: i 



lim e^^ff^o) -» 

)— t+OO 
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In the equation (|3.13|1 . for z — u, let to —> +00; for z — c, s, let <o = 0; then one 

gets 

(3.14) n u (t) = f e (t - T)LU N u (Sl)dr , 

J +00 

(3.15) n z (t) = e tL 'n z (0)+ f e {t ~ T ^ L ~~ N z {n)dr , (z = c,s). 

Jo 

The right hand side of lj3~T4) - (j3~TB) defines a map r = T u + T c + V s on B lti . First 
we show that T contracts, then T acting is an easier argument. From (|3.14(l . one 
has 

/+00 
e^ t e a{t ~ T) \\N u (Q 1 ) - N u (Q 2 )\\ e dT . 

Using H3.11JI . one gets 

/+00 
e-^e^-^W^iQi) - 7V"(0 2 )||^_irfr , 

in this article, all constants are denoted by the same C. Since H 1 ^ 1 is a Banach 
algebra 

||iV"(fix) - Ar-(0 2 )|U_! < C{\\£l x \\z + - flail* . 

Finally one gets 

/+00 
e^-^-^iWQxWt + Hnall/Je-^llfJi - ihhdr . 

Now we need to replace iV(fi) by its cut-off x(||fl|U/^)-^(£l) for some 5 > 0, where 

(3.16) X(r) = l, re [0,1]; x(r) = , re[3,+co); 

X '(r)<l, re[o,H; xeC °°(R,K). 

If both x(||fli|U/^) and x(||fl2||f/<5) are zero, then of course 

x (\\n 1 \\ e /s)N(n 1 ) - x(m 2 \\ e /s)N(n 2 ) = o . 

If one of them is nonzero, without loss of generality, say x(||fli|U/<5) 7^ 0, then 
||X(|| WWflx) - x(\\Sh\\t/S)N(Sh)\\t-i 

< ||[x(l|fii||//<*) - xiW^h/swi^We-i + ||x(||n a |U/a)[tf(ni) - N(Sh)]\\t-i 

< \x(\\^i\U/5) - x(l|fl2|k/<5)IIWi)IU-i + Ixdl^lU/^IH^Ox) - iV(Q 2 )|U_ a 

< hsii - n2||<||JV(fii)||<-i + |x(l|fi2lk/<5)|||AW - Ar(o 2 )|U_! . 



If x(||fl2lU/<5) = 0, then the second part disappears; otherwise, ||£1 2 ||^ < 3<5. In any 
case, one gets 

\\ x m 1 \\ e /6)N(n 1 ) - x{\\toa\U/S)N(Sh)\\t-i < CS\\Sl! -ih\\t. 
Thus, one has 

/+00 
e ( a _ 7 )(t_r) dT 

a — 7 
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Thus 

(3.17) supe-^llr"^! - T u n 2 \\e < C IIOi - n 3 |L t . 

For z = c in (|3.15() . one has 

e- 7t ||r c ni - r c n 2 \\ e < c [ e-^e^^WN^n^ - N c (n 2 )\\ e dT , 



here notice that fi^ and f2§ have the same initial condition fif(0) = The 
same argument as above leads to 

r V(7-e)(t-T) dT 



e-T t ||r c ni-r c n 2 || < < cs^ -n a || 7 ,/ / 

Jo 



< C HOi-fialL^. 

7- e 

Thus 

(3.i8) supe- 7 *||r c fii-r c fi 2 |b < c |ini-n a |Li . 

t>o 7-e 
For z = s in l|3.15|l . one has 

e-^||r 4 ni - T s Q 2 \\ e < C [ e-^[l3{t - T^-^e'^-^WN 3 ^) - N'(Sl2)\\t-idT 

Jo 

The same argument as before leads to 



_S_ 

v? 1 



< C— ||n x — Qall^ 



Thus 



6 



(3.19) supe-^iir-fii - r s n 2 \u < c-=\\Q! - n 2 || 7 , 

t>o VP 

By choosing 

(3.20) 5 = n 
(|3~T7|l - !j3~T5|) imply that 



m |ia(a-7),gC(7-e),^V^} , 



Thus r contracts. The claim that T acts, i.e. T : B 7t g — > -B 7j ^; follows similarly, and 
is an easier argument. By the contraction mapping theorem, T has a fixed point 
in -B 7 /, which satisfies (|3.14|) - (I3.15|) . of course, with N replaced by its cut-off. 
Equations H3.14fl - H3.15fl define a map which maps fl z (0) (z = c, s) into 

r° 

Q u (0) = / e~ TL N u {VL)dr . 

J+oo 

This map defines the center-stable manifold. Due to the cut-off, only in the 8 neigh- 
borhood of Vt* , this center-stable manifold corresponds to the NS flow. Smoothness 
can be proved through the standard argument . This proves the theorem. □ 
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To study the zero-viscosity limits of invariant manifolds, we need to make an 
assumption. 

• Assumption 1. As v — ► + , Re{A"} is of order 0(1), and Re{Af} is 
0(y), and the constant C in (f3^ - lj3~TT)l is 0(1). 

Theorem 3.6. Under Assumption 1, as v — > + , in an order Oi^Jv) neighbor- 
hood of the fixed point 11* in the Sobolev space H e (T d ) (I > 3, d — 2, 3), there exist 
an unstable manifold and a center-stable manifold; and in its order 0(v) neigh- 
borhood, there exist a stable manifold, a center manifold, and a center-unstable 
manifold. 

PROOF. Under the Assumption 1, the S in (|3.20() is of order 0(y/U) since f3 is 
0(v). Thus we have the claim for the center-stable manifold. By keeping track of 
similar estimates as in the proof of Theorem 13. 51 one can quickly see the rest of the 
claims. □ 

Due to Assumption 1, zero- viscosity limits of invariant manifolds become a challeng- 
ing problem. The challenge is further amplified by the elusiveness of the invariant 
manifolds for Euler equations as discussed in the Introduction. As a test, we studied 
a simpler problem in |12| . which reveals some unique features of the zero- viscosity 
limit of the spectra. For instance, the zero-viscosity limit of a discrete spectrum 
can be a continuous spectrum which is not the spectrum of zero-viscosity. 



4. Examples 

Using Fourier series for the 2D NS l|2.1|) , 

n = J2 ^e lk - x , / = E ^ k ' x > 

fcez 2 /{o} fcez 2 /{o} 
where uj-k — uJk and = fk, one gets the kinetic form of 2D NS 



u k = ^2 A(q, r) u q u r + v[- 

k—q-\-r 



\k\ 2 u k 



fk] 



where 



where |q| 



Afar) - \ [I? 



1x 

<12 



ri 
r-i 



q\ for q = (qi,q2) T , similarly for r. Linearize the 2D NS at the 
steady state given by a single mode, 

(4.1) n* = Fe ip - X + f e~ ivx 

where / = |p| 2 f2», one gets 



k+np 



= A(p, k + (n - l)p) T 



k-\-(n— l)p 



(4-2) + A(-p,k+(n + l)p) f u, k+{n+1)p -v\k + np\ 2 u; k+np 

where its right hand side defines the linear NS operator L. Let ^j, +np 



r = |r| e % a = i|r| 

then 
(4.3) 



Pi 

P2 



ki 

k 



, Pn = \p\ 2 - \k + np\ 2 , and z n = p n e~ 



k-\-np ' 



a n z n + Z n -1 - Z n +l = , 
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where a n = (ap n )^ 1 (X + v\k + np\ 2 ). Let w n = z n /z n -i, one gets (16) 
(4.4) a n + — = w n+1 , 



which leads to the continued fraction, 

(4.5) u>W = a„_i 4 

Rewriting (|4.4(l as, 

(4.6) w n = 



1 



a«-2 + 



1 



a n -3+" 



-a„ + w n+ i 
which leads to the other continued fraction 

1 



(4.7) w 



(2) - _. 



a„ 



a n+1 + -- 



Along the spirit of the finite difference of a second order ordinary differential equa- 
tion, the difference equation (|4.3|) should have two linearly independent solutions. 
Rigorous theory has been well developed. This is the so-called Poincare-Perron 
system. For details, see [H] HI [HI [HI [HI 0] E] H] El [HI When v > 0, 
a n — > an 2 as |n| — ► +oo. Then there are a growing solution z n /z n _i — > an 2 and 
a decaying solution z n jz n -\ — > a^ 1 n^ 2 , as n — > +oo; and vice versa as n — > — oo. 
The intuition on this is clear from 14. 4|) and (|4.6|) . The eigenvalues are then deter- 
mined by matching the two decaying solutions given by the two continued fractions 
(|4~5|) and (P) . When = 0, a„ -> a = a|p| 2 A as |n| +oo. When Re{a} ^ 0, 
or Re{a} = (|a| > 2), there are a growing solution z n /z n -x — > w+, (\w+\ > f), 
and a decaying solution z n /z n -\ — > (|w-| < f), as n — > +oo; and vice versa as 
n — * — oo. Here w± solve the characteristic equation [5] 

w 2 — o!c — f = . 

Again the eigenvalues are determined by matching the two decaying solutions given 
by the two continued fractions (|4.5(l and 1(4. 7JI . Finally, Re{a} = (|a| < 2) 
corresponds to a continuous spectrum [2] . The proof of can be generalized to H e 
setting for any natural number I, For more general result, see |25j . In summary, 
when v > or v = [Re{a} ^ 0, or Re{a} = (\a\ > 2)], the eigenvalues are 
determined by 

(4.8) a ■ 



+ „ °1 



a-2H — a2H — 

which is obtained by matching wlP (|4.5|l and (|4.7|l at n = 1, i.e. — w^K 
The corresponding eigenfunctions belong to H e for any natural number I. 



4.1. Example 1. Our first prime example is the shear fi* = cos £2 defined on 
the rectangular periodic domain [0, 2n/a] x [0, 2%] where 1/2 < a < 1. As shown 
below, this example has only one unstable eigenvalue. Complete information on 
this eigenvalue can be obtained. Thereby, existence of an unstable manifold and a 
stable manifold can be confirmed. This example is motivated by examples 2 and 3 
studied later. 
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In this rectangular domain case, Q has the Fourier expansion 

Q = ^2 LU k e l[aklXl+k2X2] . 
feez 2 /{o} 

We have equation l|4.3|l with 

4 (afci) 2 + (fc 2 +n) 2 



-{A + ^[(afc 1 ) 2 + (fc 2 +n) 2 ]} 



aki (aki) 2 + (Zc 2 + n) 2 

Theorem 4.1. The spectra of the 2D linear Euler operator L have the following 
properties. 

(1) (afci) 2 + (fc 2 + n) 2 > 1, Vn G Z. When v —> 0, there is no eigenvalue of 
non-negative real part. When v = 0, i/ie entire spectrum is the continuous 
spectrum 

.a\ki\ .a\ki\ 



2 ' 2 

(2) ki = 0, &2 = 1. TTie spectrum consists of the eigenvalues 

A = -vn 2 , ne Z/{0} . 

TTie eigenfunctions are the Fourier modes 

u np e mX2 + c.c. , Va) np € C , ne Z/{0} . 

As v — ► + , ifte eigenvalues are dense on the negative half of the real axis 
(— oo, 0]. Setting v — 0, t/ie cmZy eigenvalue is A = of infinite multiplicity 
with the same eigenfunctions as above. 

(3) /ci = — 1, fc 2 = 0. (aj. f > 0. For any a G (0.5,0.95), there is a unique 
v*{ot), 



, An . V32 - 3a 6 - 17a 4 - 16a 2 , . 1 1 

( 4 ' 9 ) Af 2 , ^/ 2 , A\ < < 



4(a 2 + l)(a 2 + 4) * v ' 2(a 2 + 1) V 2 ' 

where the term under the square root on the left is positive for a € 
(0.5,0.95), and f/ie Ze/f term is always less than the right term. When 
v > v* (a) , iftere is no eigenvalue of non-negative real part. When v = 
(a), A = is an eigenvalue, and all the rest eigenvalues have nega- 
tive real parts. When v < v m (a) } there is a unique positive eigenvalue 
\(v) > 0, and all the rest eigenvalues have negative real parts. v~ 1 \(h l ) is 
a strictly monotonically decreasing function of v. When a £ (0.5, 0.8469) , 
we have the estimate 



'a 2 (l-a 2 ) a 4 (a 2 + 3) , 2 „ , 

-8RTlT' l6(a 2 + l)(a 2 + 4) -^ + D< A M 



a 2 (l-a 2 ) , 

where the term under the square root on the left is positive for a £ 
(0.5,0.8469). 



' a 2 (l-a 2 ) _ a 4 (a 2 + 3) / a 2 (l-a 2 ) 
8(a 2 + 1) 16(a 2 + l)(a 2 + 4) " W " V 8(a 2 + 1) 
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In particular, as v — > + , \(v) = 0(1). 

(b). v — 0. When a € (0.5,0.8469), we have only two eigenvalues Ao 
and — Aq 7 where Aq is positive, 



l a 2 (l-a 2 ) a 4 (a 2 +3) ^ < /a 2 (l 



V 8(a 2 + l) 16(a 2 + l)(a 2 + 4) u y 8(a 2 + 1) ' 

TTie rest of the spectrum is a continuous spectrum [—ia/2, ia/2]. 
(c). For any fixed a € (0.5, 0.8469), 

(4.10) lim \(v) = A . 

v— >0+ 

(4) Finally, when v = 0, £/ie union of all the above pieces of continuous spectra 
is the imaginary axis iM.. 

From Theorems 13.51 and 14 . 1 1 we have the corollary. 

Corollary 4.2. For any a S (0.5,0.95), and v S (0,M*(a)) w/iere f*(a) > 
satisfies \4-^ , in a neighborhood of Vl* in the Sobolev space H e (T 2 ) (I > 3), there 
are an l-dimensional C°° unstable manifold and an 1-codimensional C°° stable 
manifold. 

Remark 4.3. In the Theorem 14. II I4.10|l verifies part of Assumption 1. That 
is, Re{A"} is of order O(l) as v — > + . Case 2 of Theorem l4.1l indicatcs that Re{Af} 
is at least 0(v) as v — > + . Case 3 indicates that the constant C in l|3.7|l and (|3.11() 
should be 0(1) as v -> 0+. 

Proof. We will give the proof case by case. 

(1) The case (a), v > 0. If Re{A} > 0, then all the Re{a n }'s are nonzero 
and have the same sign. The real parts of the right and left hand sides 
of (|4.8|l are nonzero but of different signs. Thus there is no eigenvalue 
of non-negative real part. The case (b). v = 0. If Re{A} ^ 0, then all 
the Re{a ra }'s are nonzero and have the same sign. Again (|4.8(l can not be 
satisfied. If RejAl = 0, let a = lim„^co a n = — then 

L J aki 

_ (aki) 2 + (h + n) 2 - 1 

az n — a -~ — (Zn+l - Zri-lj j 

(afci) 2 + (fe 2 +n) 2 

where (aki) 2 + (k% + n) 2 > 1, Vu e Z. By using £2 norm of {z ra }n6Z, 
one sees that possible eigenvalues have to satisfy \a\ < 2. As mentioned 
before, Re{A} = and \a\ < 2 correspond to a continuous spectrum |9], 
which is the interval 

a\ki\ a\ki\ 
-i , 1 



(4.11) 

(2) In this case, one has 

[A + v(n + l) 2 ]d)(„ + i)p = . 

The claims follow immediately. 

(3) In this case, 

4 a 2 + n 2 o 2 

a n = , 5 A + v(a +n 

a a z + n l — 1 
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Thus a_„ = On. Equation (|4.8(l is reduced to 
(4.12) -a /2 ' 



ai 



The first a few a„'s are 

&o = — ^ [A + va \ 

ai = ~ ~ „.2 ^ [A + ^(q 2 + 1)] 
fl2 = ' ^ [A + ^(a 2 +4)] 



Let 



4 a 2 


a a 2 — 


1 


4a 2 + 


1 


a a 2 




4a 2 + 


4 


a a 2 + 


3 


1 




j. 1 



/(A) = ; — , g(X) = -a /2 

a-i 



a-2-\ — 

(a), v > 0. First we show that if A (Re{A} > 0) is a solution of (|4.12[) . 
then A must be real. Assume Im{A} > 0, then 

arg(-ao) > arg(ai) > arg(a 2 ) > • • • > . 

Thus 

| arg(/(A))| < arg(oi) . 

But 

arg(g(A)) = arg(-a ) > | arg(/(A))| , 

a contradiction. Similarly for the case Im{A} < 0, thus A (Re{A} > 0) 
is real. Next we show that there is a vo > such that for every v < vo, 
there is a unique eigenvalue A > 0. Since a n > 0, Vn > 1, we have 

(4.13) — ^ < /(A) < - . 

Thus when A is large enough 

(4.14) /(A) < 5 (A) . 
From (|4~T3|) . 

1 



c*3 + 4y(a 2 +4) 2 

We know that 

2z/a 3 
1 — or 

We want to choose ^ such that 

1 2isa 3 

> 



4t/(q 2 + l) 2 aQ 2 +3) - \ _ q,2 

h 4i^(a 2 +4) 2 



Thus 



1C , ^ , . V32 - 3a 6 - 17a 4 - 16a 2 f1 .„ . 

4.15 i/ < i/„(a) = 1/2 < a < a 

4(a^ + l)(ar + 4) 

where 0.95 < ao < 1 and «o satisfies 

32 - 3a 6 - 17a 4 - 16a 2 = . 
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For example, 

i/ (0.5) « 0.244 , i/ Q (0.95) » 0.0329 . 
For each fixed z/, v < uo(a), we have 

(4.16) /(0) > 5 (0) . 

From l|4.14|) and (|4.16|) . we see that there is a A > such that (|4.12|) is 
true. Next, we want to show that this eigenvalue is unique. Notice that 

2<T 

(4.17) (\ + va 2 )- l g(\) = ~ 



1-a 2 ' 

(4.18) (A + ^ 2 )- 1 /(A) : ' 



(A + vc?)a x + (XW) 1 - 1 ^^-- ' 

Since (X+va 2 )a2n+i (n > 0) is a strictly monotonically increasing function 
of A for A > 0, and (A + vo?)~ 1 a2n (n > 1) is a strictly monotonically 
decreasing function of A for A > 0, we see that (X + va 2 )^ 1 /(A) is a strictly 
monotonically decreasing function of A for A > 0. Thus the eigenvalue 
which satisfies 

(A + ua 2 )~ 1 f(X) = (A + ua 2 )^ 1 g(X) 

is unique. Similarly, va 2 a,2 n +i {n > 0) is a strictly monotonically increas- 
ing function of v for v > 0, and (va 2 )~ 1 a2 n ( n > 1) is a constant function 
of v. Then (va 2 )^ 1 f(0) is a strictly monotonically decreasing function of 
v for v > 0. We know from above that when v = i/o(a), 

2a 

(i/a 2 ) _1 /(0) > (^a 2 ) _1 3(0) = ? , constant in i/ . 

1 — a 2 

From (|4.13|) . we have 

{va 2 )- 1 !^) < 



(va 2 ) ai A{a 2 + l) 2 v 2 ' 

thus when v > is large enough 

(i/aTVW < (m 2 )- 1 ^) . 

Therefore, there is a unique i/*(a) > 1/0(0;), such that 

(4.19) (**o: 2 r7(0) = (v^r'giO) . 

Thus we have shown the following claims: When v > 14(a), there is 
no eigenvalue of non-negative real part. When v = i/*(a), A = is an 
eigenvalue, and all the rest eigenvalues have negative real parts. When 
v < i/*(a), there is a unique positive eigenvalue X(v) > 0, and all the rest 
eigenvalues have negative real parts. An estimate for can be obtained 
from (|4"T3|) . 

1 2a 1 



< 7; < 



4(« 2 +l) 2 ,.2 I a 3 (a 2 +3) 1 _ a 2 ^ 4(a 2 +l) 2 ,, 2 
a ^* 4(a 2 +4) 2 a ^* 

We have 

1 



z/o(a) < i/* (a) < 



2(a 2 + 1) 
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(4.20) 



(4.21) 



(4.22) 



where Vq{q) is given by (|4.15(l . For example, 

0.244029 < z/*(0.5) < 0.244949 , 0.0329 < ^(0.95) < 0.058 . 

Next we want to show that the unique eigenvalue A(i^) has the property 
that ^ _1 A(^) is a strictly monotonically decreasing function of v. Notice 
that 

>,2 _i_ (O^, _L 1)2 



(A + VOL )d2 



n+1 



(2n 



a a 
X 
v 



(X + va 2 ) 1 a 2n = -- 



2 + (2n + l) 2 - 1 
+ [a 2 + (2 ? i + l) 2 ; 

a 2 + (2n) 2 



A 



(n > 0) , 



(2n) 



1 



1 



(2n) 2 
- +a 2 



(n > 1) 



Assume that v~ X(u) is not a strictly monotonically decreasing function 
of v, then there is an interval in which v~ 1 X(u) is a strictly monotonically 
increasing or a constant function of v. In that interval, \X{y) + vot 2 ]d2n+i 
(n > 0) is strictly monotonically increasing, and [X(v) + va 2 ]~ 1 a,2n (n > 1) 
is monotonically decreasing. Thus \X(v) + va 2 ] -1 f (X(v)) is a strictly 
monotonically decreasing function of v. On the other hand, 

2a 



[\{v) + va^g{\{v)) = 



1 - a 2 



is a constant function of v. This contradiction shows that v 1 X(v) is a 
strictly monotonically decreasing function of v. Next we want to show 
that X{v) = 0(1) as v -> 0+. From (|4~T2j) and we have 

1 

4Q 2 +1)[A+Ka 2 + 1)] q(q 2 +3) 

h 4(Q 2 +4)[A+iy(a 2 +4)] 



< 



Let -> H 



2a 



1 — a 
, we have 



(A + va 2 ) < 



4(a 2 + l)[A + ^(a 2 + l)] 



'a 2 (l-a 2 ) 



a 4 (a 2 +3) 



1) 



16(a 2 + l)(a 2 +4) 



< A < 



/a 2 (l-a 2 ) 
8(a 2 + 1) 



For the term under the square root on the left to be positive, we need 

I w 0.8469 . 



a < ai 



Thus when a G (1/2, ai), 

A(f) = 0(1) , as v -> H 
Also from (|4.20|) . we have the estimate 



'a 2 (l 



a 4 (a 2 +3) 



8(a 2 + l) 16(a 2 + l)(a 2 +4) 



- z^(a 2 + 1) < \{v) 
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(b). u = 0. We have 



of 



(4.24) 



-A . 



ft 9 i 9 1 

As \n\ — * oo, a n — ► a = -A. As before, Re{A} = and |A| < § cor " 
responds to a continuous spectrum. Next we show that outside the disc 
|A| < §, there is no eigenvalue. Outside the disc, \a\ > 2, 

\a a \ > 2^-^ > I , force (1/2,1) , 
z 3 



1 



a- 



1 



|a 1 |>2ll + -l>4 ) for a 6(1/2,1) 
and 

\a„\ > 2 , Vn > 2 . 
By a result of continued fraction |15| . 

1 



a 2 



a 3 + - 



< 1 



Thus 



while 



\gW\ > 



l/(A)| < 



»3 + " 



< 



This contradiction proves the claim. Thus the possible eigenvalues should 
lie in the region Re{A} ^ and |A| < Since the eigenvalues are 
symmetric with respect to the imaginary axis [5] , we only need to consider 
the case Re{A} > 0. First we show that the possible eigenvalue must be 
real. Assume that Im{A} > 0, then 

arg(-a ) = arg(a x ) = arg(a 2 ) = • • • > . 

Thus 



Then 




arg ai 



Therefore 



arg(/(A)) jt arg(oi) = arg(-a ) = arg(.g(A)) . 

This is a contradiction. The case Im{A} < is similar, and the claim is 
proved. When A > 0, 

■3— < A-\f(A) < -L, 



Aai 
1 



Aai 



4(a 2 + i; 
^3 



A 2 



a(a 2 +3) 
4(a 2 +4) 



< A-7(A) < 



1 



4(« 2 + l) A2 
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Thus 

lim A- 1 /(A) = 0, lim A-V(A)> 4( " 2 + ^ 

A^+co A^0+ a{a z + 3) 

Notice that 



1 

We want to choose a such that 



2a _ 4(a 2 + 4) 



1 — a 2 a(a 2 + 3) 

This condition is equivalent to the fact that the term under the square 
root on the left of Q4.21I) being positive. Thus we have 



a < at = \ J w 0.8469 , 

V V 12 2 

which is the same with (|4.22|l . Thus when a G (1/2, ai), there is a positive 
eigenvalue Ao, 

Aq V(Ao) = Aq ^(Ao) = y_~~2 • 

Since \ci2n+i {n > 0) is a strictly monotonically increasing function of A 
for A > 0, and \~ 1 a,2 n (n > 1) is a constant function of A, we see that 
A _1 /(A) is a strictly monotonically decreasing function of A for A > 0. 
Thus the positive eigenvalue Ao is unique. From (|4.24|l . we have 
1 2a 1 

< n < 



4(a 2 + l) x2 | a(a 2 +3) 1 - a 2 4(q 3 +l) x 2 ' 
^? A ' 4(a 2 +4) a 3 A 



i.e. 



a 2 (l-a 2 ) a*(a 2 +3) a 2 (l-a 2 ) 

[ ' V 8(a 2 + 1) 16(a 2 + l)(a 2 + 4) \j 8(a 2 + 1) ' 

which is similar to 14. 2111 . Next we want to show that 

(4.26) lim \{v) = A . 

Let F(v, A) = /(A) — g(A), and let Fn(i/,X) be the 7V-th truncation of 
F(p,X), 

F N (v,\) = ^+ ■ 1 - . 

2 "1 + a2 +- + sL 

Notice that 

F(^,A(^)=0, ^ G (0,J/*(a)) , a G (0.5, 0.8469) ; 

F(0,A ) = 0, a G (0.5, 0.8469) ; 

where A(i/) and Ao have the estimates (I4.23|) and l|4.25|l . Thus for any 
fixed a G (0.5, 0.8469), (v, lies in a compact set [0, S] x [ci, c 2 ] where 

S, ci and C2 are positive constants. Assume that H4.26J) is not true, then 
there is a sequence (yj, such that 

lim \{vj) = Ao , 

Vi->Q+ 
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where 

A ^ A and A G \c\ , c 2 ] . 
By a result of continued fraction 

lim Fn{v, A) = F(v, A) 

uniformly for z/ € [0,5] and A e [ci,C2]. Thus for any e > 0, there is an 
integer Nq such that for any N > iVo, 

[Fjv^^A^))-^^^))! <e, 

i.e. 

For any fixed N, let j — > +oo, we have 
|i ? V(0,Ao)|<e. 

Let TV — > +oo, we get 

1^(0, A )| <e • 
Since e is arbitrarily small, we have 

F(0,A o ) = . 

By the uniqueness of the eigenvalue, Ao = Ao. This contradiction shows 
that (|4~2ll is true. 
(4) When v — 0, the union of all the continuous spatra 

alfcil a\ki\ 
-i , i 



□ 

4.2. Example 2. Our second prime example is = cos(.ti + X2) which 
corresponds to p — (1, 1) T and r = 1/2 in 14.1JI . Here 



(fci +n) 2 + (k 2 + n) 2 



(k 2 - fci) (fci + n) 2 + (k 2 + n) 2 - 2 



(4.27) 



x < A + v 



(fci 



(fc 2 + n) 2 ]} 



This example has two possible unstable eigenvalues. For shears, the number of 
unstable eigenvalues of 2D linear Euler is bounded by the number of lattices points 
inside the disc of radius \p\ [8], which is even for square lattice. In |13| |14j . Vincent 
Liu studied another shear which has three possible unstable eigenvalues. Liu did a 
detailed calculation on the eigenvalues, which will be discussed in next subsection. 
Here we will do an even more detailed calculation. 



Theorem 4.4. The spectra of the 2D linear Euler operator L have the following 
properties. 
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(1) The set {k + np} ne z has no intersection with the disc of radius \p\ = \pl. 
When v > 0, there is no eigenvalue of non-negative real part. When v = 0, 
the entire spectrum is the continuous spectrum 

-i-\k 2 - ki\, i^\h - fci| ■ 

(2) k = p = (1, 1) T . The spectrum consists of the eigenvalues 

X = -v2n 2 , n S Z/{0} . 

The eigenf unctions are the Fourier modes 

Cj np e mp ' x + c.c. , yCj np E C , n € Z/{0} . 

As v — > + , i/ie eigenvalues are dense on the negative half of the real axis 
(— oo, 0]. Setting v — 0, the only eigenvalue is X — of infinite multiplicity 
with the same eigenf unctions as above. 

(3) k = (-1,1) T . When v > 0, there is no eigenvalue of non-negative 
real part. When v — 0, the entire spectrum is the continuous spectrum 
[— iijig]. j4 special eigenvalue is X = —2v ( when v = 0, i/iis eigenvalue 
X = is embedded in the continuous spectrum). 

(4) k = (0, 1) T . When v > 0, in the half plane Re{X} > —v, there is a unique 
pair of eigenvalues X and X such that 

-v< Re{X} <yi + (8vy-2v, 

I(-i/I)<^k|( 1+ ^)- 

When v = 0, [— ij, ij\ is a continuous spectrum. If there is an eigenvalue 
of positive real part, then there is a quadruplet (X, X —A — X) where 

(5) Finally, when v = 0, the union of all the above pieces of continuous spectra 
is the imaginary axis iM.. 

Remark 4.5. For Case 4, numerical computation indicates that when v > 0, 
the real part of the eigenvalue Re{A(^)} > c > 0, where c is independent of v. When 
v = 0, numerical computation indicates that the eigenvalue A(0) indeed exists and 
Re{A(0)} > c > (the same c as above). Numerical computation also indicates 
that as v — > + , A(^) — > A(0) which can be proved given the above facts. Cases 2 
and 3 indicate that Re{Af } in Assumption 1 is at least 0(y) as v — * + . 

Proof. We will give the proof case by case. 
(1) The case (a), v > 0. If Re{A} > 0, then all the Re{a„}'s are nonzero 
and have the same sign. The real parts of the right and left hand sides 
of (|4.8|l are nonzero but of different signs. Thus there is no eigenvalue 
of non-negative real part. The case (b). v = 0. If Re{A} ^ 0, then all 
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the Re{a n }'s are nonzero and have the same sign. Again l|4.8(l can not be 
satisfied. If Re{A} = 0, let a — lim„^co a n = ; 8 ; A, then 

\k + np\ 2 -2 



np\ 2 



-[z n +l ~ ^n-l) ! 



where < \k + np\ 2 — 2 < \k + np\ 2 , Vn G Z. By using £2 norm of {z n } n ^x, 
one sees that possible eigenvalues have to satisfy \a\ < 2. As mentioned 
before, Re{A} = and \a\ < 2 correspond to a continuous spectrum [9], 
which is the interval 



(4.28) 

(2) In this case, one has 



ij\k 2 - h\, - fci| 



[A + v2(n + 1) ]iD( n+ i)p = . 

The claims follow immediately. 
(3) In this case, denote &% +np simply by w n , one has 

n , 2l1 „ l (n + l) 2 l (n-l) 2 

[A + 2,(n 2 + l)H = 4 (B Vi )a + i ^ - I (n -l)2 + l ^ 



This system decouples at n = 0. A few equations around n = are 

9 

40 C 



. . 19 

(A + 10^)oj_2 = -u-i - 77:^-3 , 



(4.29) (A + 4i/)w_i = --w_ 2 , 

5 

1 1 

(A + 2^)a> = g^l - gW-i , 

(4.30) (A + 4i/)wi = , 

9 1 
(A + 10^)^2 = 77:^3 — a UJl ' 



40 



Notice that 



1 



4 1 + — [\ + 2v{n z + 1) 



From (|4.30|) . one gets 

(4.31) w 2 = z 2 /z! = 8(A + Av) 
From (|4.5|l and (|4.31|) . one has 

(4.32) 8(A + 4i/) ' 



a-2 



(a). 1/ > 0. If Re{A} > 0, then Rc{a„} > for n > 2. Thus the real 
part of the right hand side of 14.32J1 is negative, while the real part of 
its left hand side is positive. Thus there is no eigenvalue of non-negative 
real part. (b). v = 0. If Re{A} > 0, then Re{a„} (n > 2) has the same 
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fixed sign with A, again l|4.32|) leads to a contradiction. If Rc{A} = 0, let 

a = lim n _^oo a n = 4A, then 



1 



-(z n+ x - z n _i) , n > 2 



n 2 + 1 

By using £2 norm of {z n } n >i, one sees that possible eigenvalues have to 
satisfy \a\ < 2. As mentioned before, Re{A} = and \a\ < 2 correspond 

i± ail 

'2 ' 2J 



to a continuous spectrum which is the interval [— Similarly 



from (|4.29l) . one gets 

1 



(4.33) w~i=Z-i/z-2 
From (|4.5|l and (|4.33() . one has 



8(A + Av) 



1 1 

a- 2 



8(A + Av) a_ 3 



a_ 4 + --- 

Notice that a_„ = a„, one gets the same conclusion as above. Finally, 
by choosing u> n = if n 7^ 0, one gets the eigenvalue A = —2v with the 
eigenfunction 



ik-x 

(4) In this case, 



u~ k e + c.c. , Muif, e 



+ (n + \y — 2 
Thus a_( Jl+1 ) = a n . Equation (|4.8|) is reduced to 

(4.34) a Q + — p- = ±* . 

The first a few a n 's are 

a = 8-(-l)-(A + i/) , 

ax = 8 ♦ (|) - (A + 5i/) , 

a 2 - 8 -(H). ( A + 13!/) , 
25 

«3 = 8 ■ (-) ■ (A + 25*) . 

(a), v > 0. Consider the region Re{A} > — 1/, in which 

\a n \ > 8\X + v\ , Vn > . 

By using ^2 norm of {z n }nez in (|4.3(l . one sees that possible eigenvalues 
in this region have to satisfy 

\\ + v\<\. 

Next in the region 

1 



V = { A 



Re{A} > -v, |A + v\ < 4 
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we will use Rouche's theorem to track the eigenvalues. Let 

1 



/(A) 



ao + 



ai 
1 



a2+- 



g(X) = a H hi. 

ai 

To apply the Rouche's theorem, one needs to show that 

|/(A)- fl (A)|<|/(A)| + |fl(A)|, XedV. 
If this is not true, then there is a A € &D and a 6 (0 < S < oo) such that 
/(A) = -5g(X) , 



a + i + 



l 

1+5 



1 



«1 



1 + <5 ai 







On the part of the boundary dT>: Re{A} = — v, 

Re{a } = , Re{a„} > , n > 1 . 

By taking the real part of (|4.37() . one sees that (|4.37() can not be satisfied. 
On the other part the boundary dT>: |A + v\ = j, 



\a \=2, \ ai \> 



10 



M > 



26 
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\a n \ > 2 , Vn > 3 . 
By a result of continued fraction |15) . 

1 



«3 



< 1 



Thus 



l 

1+5 



a i 



02 + " 



< 



< 



< 



ax 



0.2 + - 



ax 



10 1 
"3 W- 



>3 + ^ 

15 
" 39 ' 



The rest of (|4.37|l has the estimate: 

S 1 



a + i + 
> 2- 1- 



1 



1 + 5 ai 
5 3 
l + JlO " 10 



1 + 5 bi 



> 
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Therefore, Ij4.37|) can not be satisfied. By Rouche's theorem, inside T>, 
/(A) and g (A) have the same number of roots. g(X) has two roots 



A 



(4.38) 



and one of which is in D (at least when v is small enough). Thus /(A) 
has one root in T>. From (|4.35[l . this root satisfies 

1 



Re <y a\ + 



a + i 



<0 , 



which leads to 



-v < Re{A} < 



20 



(81/) 



2v 



< Im{A} <^1 + Vsj • 
(b). v = 0. Then 

\a„\ > 8|A| , Vn > . 

Again by using i 2 norm of {z n } ne z in l|4.3|l. one sees that possible eigen- 
values in this region have to satisfy 

1 



IAI < 



In fact, if there is an eigenvalue with positive real part, then 
which leads to 

o< ma} <iyf , 



is true, 



1 



< lm{A} < if 1 



As proved in [S], such eigenvalues (if exist) come in quadruplet (A, A —A 
— A). As in Case 3 above, [— ij, ij] is a continuous spectrum 
(5) When v = 0, the union of all the continuous spectra 

-i-\k 2 - fei|, - 



is 



□ 



4.3. Example 3. An example studied in details by Vincent Liu |13) |14| is 
= — ^ cos(2x2) which corresponds to p = (0, 2) T and T — in (|4.1() . Here 



kl 



(fe 2 + 2n) 2 



k\ + (k 2 + 2n) 2 - 4 



{A + v \kj + (fc 2 + 2n) 2 ] } 



Theorem 4.6. The spectra of the 2D linear Euler operator L have the following 
properties. 
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(1) The set {k + np} n ez has no intersection with the disc of radius \p\ = 2. 
When v > 0, there is no eigenvalue of non-negative real part. When v = 0, 
the entire spectrum is the continuous spectrum 

■ M . \k\ 



2-V/27I-' 2s/2ir 

(2) k = (0, 1) T or (0, 2) T . T/ie spectrum consists of the eigenvalues 

A = + (fc 2 + 2n) 2 ] , neZ, 

w/iere /or fc = (0, 2) T , n ^ — 1. TTie eigenf unctions are the Fourier modes 

*k +np e iik+np>x + ** , Vu> fe+np € C , » e Z . 

is i/^ + , £/ie eigenvalues are dense on the negative half of the real axis 
(—oo, 0]. Setting v = 0, the only eigenvalue is A = of infinite multiplicity 
with the same eigenf unctions as above. 

(3) fc = (2,0) T . Mien z/ > 0, t/iere is no eigenvalue of non-negative real part. 
When v = 0, i/ie entire spectrum is the continuous spectrum . 
A special eigenvalue is A = —4^ (when v = ; t/iis eigenvalue A = is 
embedded in the continuous spectrum). 

(4) fc = (1,0) T . T/iere is a unique v*, 

1 /89 1 /— 

IOtt V 34 IOtt 
When v > v*, there is no eigenvalue of non-negative real part. When 
v = v*, A = is an eigenvalue, and all the rest eigenvalues have negative 
real parts. When v < v*, there is a unique positive eigenvalue X{v) > 0, 
and all the rest eigenvalues have negative real parts. v \{y) is a strictly 
monotonically decreasing function of v. \{v) has the estimate 

1 - w x 1 A3~ 

— \ —— - hv < X(v) < —\— — v . 

7T V 680 y ' 7T V 20 

In particular, as v — > + , X(y) = 0(1). When v = 0, we have only two 
eigenvalues Xq and —Xq, where Xq is positive 



1/89 1/3 
^V 680 <A0< ^V^0 • 

The rest of the spectrum is a continuous spectrum More- 
over, 

lim X(v ) = A . 

(5) k = (1,1) T . When v > 0, in the half plane Re{X} > —2v, there is a 
unique pair of eigenvalues X and X such that 

-2v < Re{X} < - 

IT 

V5-V3 V5 + V3 

= — < Im\X\ < = — 

2VT0tt 2VT0tt 
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When v = 0, \—i — \^,i — 7=-] is a continuous spectrum. If there is an 

' L 2V2tt' 2V2ir' e J _ 

eigenvalue of positive real part, then there is a quadruplet (X, X —A — X) 
where 

< Re{X} < 1 ; ' ! 

— — < Im{X\ < = — . 

(6) Finally, when v = 0, the union of all the above pieces of continuous spectra 
is the imaginary axis iM.. 

Remark 4.7. Except for Case 4, the rest of the proof is similar to that of 
Theorem 14.41 Overall, Liu |13j |14| did not realize the continuous spectrum. For 
Case 4, Liu |13| |14| gave an elegant proof. For Case 5, Liu |13j |14| claimed 
more than what were actually proved. In |14| . the sketch after line 9, pp.472 can 
be realized as shown in the proof of Case 4 of Theorem 14.41 above. But Equation 
(39) on the same page does not imply Re{r]i(v)} > in Equation (5) on pp.467. 
Arguments between lines 5 and 8 on pp.483 are not solid, therefore, the existence 
of r)i in Part B of Theorem 2 on pp.467 was not proved. The sketch between lines 
22 and 25 on pp.484 is not completed, therefore, (14) in Theorem 3 on pp.468 was 
not proved. Proving these claims seems tricky. 

Remark 4.8. Cases 2 and 3 indicates that Re{Af } in Assumption 1 is at least 
0(v) as v — ► + . Case 4 shows the existence of an unstable eigenvalue. In fact, 
by the property of continued fraction, the corresponding eigenfunction Q u (y) also 
converges as v — ► + . Thus ||f2 u ||^ + i/||f2 u ||^ has a uniform bound as v — > + . 

Remark 4.9. In an effort to minimize the number of unstable modes to one and 
to have Case 4 not Case 5 in the above theorem, we need to study the rectangular 
periodic domain [0, 2n/a] x [0, 27r] where 1/2 < a < 1, and the steady state fi* = 
cos £2- This is our Example 1 studied before. 
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